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ABSTRACT 

Encoded within the morphological structure of galaxies are clues related to their for- 
mation and evolutionary history. Recent advances pertaining to the statistics of galaxy 
morphology include sophisticated measures of concentration (C), asymmetry (A), and 
dumpiness (S). In this study, these three parameters (CAS) have been applied to a 
suite of simulated galaxies and compared with observational results inferred from a 
sample of nearby galaxies. The simulations span a range of late-type systems, with 
masses between ^10^° M© and ~10^^ Mq, and employ star formation density thresh- 
olds between 0.1 cm~'^ and 100 cm~^. We have found that the simulated galaxies 
possess comparable concentrations to their observed counterparts. However, the re- 
sults of the CAS analysis revealed that the simulated galaxies are generally more 
asymmetric, and that the range of dumpiness values extends beyond the range of 
those observed. Strong correlations were obtained between the three CAS parame- 
ters and colour (B-V), consistent with observed galaxies. Furthermore, the simulated 
galaxies possess strong links between their CAS parameters and Hubble type, mostly 
in-line with their observed counterparts. 
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1 INTRODUCTION 

Classification in astronomy has led to many advances which 
have revealed important information about the Universe in 
which we live. For example, the classification of stars by their 
colours and brightnesses in the Hertzsprung-Russel (HR) 
diagram led to an understanding of stellar structure and 
evolution. The HR diagram continues to be employed as a 
means to to determine the ages and metallicites of stellar 
populations, whether they are simple stellar populations, 
as in globular clu sters, or in composi t e populations , such 
as dwarf galaxies JDolphin et al.] l2003l : iMonelli et "afl I2OI0I : 
IWilliams et al.ll201Cl ). Even when it is impossible to resolve 
stars, entire gala:xies are classified by their colours and mag- 
nitudes. However, galaxies display resolved structure and 
morphology that stars do not, and using this morphology 
can reveal much about the evolution of galaxies. 
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Early methods used for th e classifi c ation of galaxies 
were based on morphology. The iHubbld (| 19261 ) Classifica- 
tion System is a method used to categorise galaxies by their 
morphology, and although it has been the dominant mor- 
phological tool since the mid-1920s, over time it has proven 
to be deficient in several areas. A prominent disadvantage 
to classifying galaxies based on their morphological features 
alone is that it is subjective with respect to distance (or 
resolution) and inclination. Furthermore it has led to the 
grouping of a wide range of asymmetric galaxies as simply 
"irregular" . 

As both ground- and space-based imaging has im- 
proved, and we probe higher redshifts, it has become appar- 
ent that classifying galaxies as "irregular" means ignoring a 
vast a mount of morphological inform ation, de VaucouleursI 
II953) attempted to improve upon the lHubblg l|l926l l Classi- 
fication System by introducing "later types" as a sub-class of 
galaxies; however, this system still refers to local, axisym- 
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metric galaxies as reference points, and, again, was based 
solely on morphology. 

[Morgan fc Mavalll () 19571 ) noted a correlation between 
the spectra of gal axies and the d ominance of their central 
bulge component. iMorganI (119581 ) felt that a shortcoming 
of the Hubble classification scheme was that it was not a 
strong indicat or of spectral class. In order to rectify this, 
iMorganI l| 19581 ') devised a galaxy classification system based 
upon the central light concentration compared to its overall 
light distribution. This concentration classification scheme 
involved analysing the flux of galaxies to determine the de- 
gree of concentration of the central component. The purpose 
of this scheme was to study the colour-concentration corre- 
lation in an attempt to identify galactic evolutionary states. 
Due to the concentration parameter's direct relation with 
spectral class and type of stellar population, the concentra- 
tion parameter was also found to be indicative of formation 
history and properties, such as velocity dispersion, galaxy 
size, luminosity, and (more recently) central black hole mass 
IJGraham et al.ll200ll '). 

E xtending Morgan's pioneering work, IConselice et al.l 
l| 20001 1 proposed a new statistical measure of asymmetry, in 
conjunction with a revised parametrisation of Morgan's con- 
centration index. This classification system was based on the 
desire to classify galaxies quantitatively at a range of red- 
shifts. Asymmetry was defined by rotating a galaxy 180°, 
about a central axis and finding the absolute sum of the 
normalised residuals. This method measures the high fre- 
quency structure within each galaxy, whilst also considering 
the symmetry. The measure of asymmetry focuses primarily 
on morphological shape, and has a strong correlation with 
colour, congruent with Morgan's concentration parameter. 
As such, this property has a correlation with both star for- 
mation and merger history. 

IConselicd (120031 ) introduced one additional morpholog- 
ical measure sensitive to high spatial frequency dumpi- 
ness. This parameter is defined by comparing a galaxy to 
a smoothed image of itself. In doing so, the dumpiness pa- 
rameter quantifies galaxy morphology in a manner which 
reflects star forming regions and evolutionary history. This 
parameter distinguishes between varying types of galaxies 
and thus, partially, resembles Hubble's classification sys- 
tem. However, the dumpiness parameter quantifies galaxy 
morphology with respect to intrinsic characteristics and not 
morphology alone. 

IConselicd (|2003l ) combined the high spatial frequency 
dumpiness parameter with both the asymmetry index and 
the concentration index to complete what is now called the 
CAS classification system. Underlying the CAS system's em- 
pirical nature lies the fundamental physics which governs its 
individual C, A, and S, components. The CAS system high- 
lights the intrinsic nature of galaxies by considering their 
light distributions and spectral class. 

The observable properties of galaxies are a function of 
their merger histories, masses and environments, each of 
which are reflected in the CAS system. It is fitting to apply 
this new classification method to the theoretically predicted 
properties of galaxies and use it as a gauge to determine 
where our current understanding of galaxy formation and 
evolution is correct and where it is deficient. Galaxy simula- 
tions are an integral tool that can help interpolate between 




Figure 1. Left panel: NGC 4216, an edge-on disk galaxy from 
the Frei et al. (1996) sample in the SDSS r-band. Right panel: 
NGC 4254, a face-on disc from the Frei et al. (1996) sample in 
the SDSS r-band. 



known stages of formation history, and thus further our un- 
derstanding of galaxy evolution. 

Within the observational community, the CAS sys- 
tem has become one of the most widely-used mea- 
sures of galaxy morpholog y (|Hernandez-Toledo et al.l [20081 : 
iBertone fc Conselicell2009^ . Since the completion of the com- 
bined CAS system, it has been applied t o numerous galax - 
ies including the 113 galaxies observed bv lFrei et al.l (|l996l ). 
However, it has not yet been applied systematically to high- 
resolution computational galaxy simulations as a method 
of comparison between observations and simulations. In 
this work, we d e velop a sof tware package pa tterned on the 
IConselice et~aLl (|2000l ) and IConseUcd (|2003l ) CAS system, 
calibrated on an optimised training set, and applied to high 
resolution galaxy simulations. We show that the simulated 
galaxies exhibit both similarities and differences to real 
galaxies, and highlight the intrinsic physical process that 
are responsible for the differences. 

In ij2]we describe our simulated galaxy samples. In ij3]we 
discuss the three structural parameters of the CAS system. 
Concentration (C), Asymmetry (A) and dumpine ss (S). We 
then c ompare the simulations with the empirical iFrei et al.l 
(|l996l ) sample of galaxies in iQ] The results are discussed in 
^and conclusions drawn in i|6] 



2 SIMULATED GALAXY SAMPLE 

To compare with the morphological parameters inferred 
from observed galaxies, we used three samples of simulated 
galaxies. Each sample was generated using the gravitational 
N-Body -I- smoothed particle hydrodynamics (SPH) code 
GASOLINE (IWadslev et ahllJOM ). We briefly outline the sim- 
ulations here, but refer the reader to the detailed descrip- 
tions found in lBrooks et al.l (120091) . for what will he nceforth 
be called the "UW sa mple" , Governato et al.l l|2010l ). for the 
"Dwarf sample" , and IStinson et al.l (|2010l ) for the "MUGS 
(McMaster Unbiased Galaxy Simulations) sample". 

Regardless of the aforementioned 'sample' from which a 
specific simulation was drawn, each individual "zoom-style" 
simulation was run within a WM AP-3 ACDM cosmologi- 
cal framework (jSpergel et al.l 120071 ) and consisted of a cen- 
tral high resolution region (centred on the target halo) em- 
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Figure 2. Left panel: An image of gl536, a disc galaxy from 
the MUGS sample, as viewed edge-on in the SDSS r-band. Right 
panel: gl536, a disk galaxy from the MUGS sample, as viewed 
face-on in the SDSS r-band. 



Figure 4. Left panel: An image of DG3, a disc galaxy from the 
Dwarf sample, as viewed edge-on in the SDSS r-band. Right panel: 
An image of DG3, a disk galaxy from the Dwarf sample, as viewed 
face-on in the SDSS r-band. 




Figure 3. Left panel: An image of MW, a disc galaxy from the 
UW sample, as viewed edge-on in the SDSS r-band. Right panel: 
An image of MW, a disc galaxy from the UW sample, as viewed 
face-on in the SDSS r-band. 



bedded within a lower resolution cosmological volume. Star 
formation and supe rnova feedback was c omputed using the 
recipes outlined by IStinson et al.l l|2006r ). with the primary 
difference between the three samples being the adopted star 
formation density threshold (nth). For the UW sample, the 
canonical nth of 0.1 cm~^ was employed; for the MUGS 
galaxies, nth was increased to 1 cm~'^; for the Dwarf galaxy 
sample, nth = 100 cm~'^ was adopted. The latter simulations 
represent the first successful realisations of essentially bul- 
geless disc galaxies (Governato et al. 2010). The masses of 
the MUGS and UW samples are comparable to that of the 
Milky Way (~10'^^ Mq) and, therefore, also the representa- 
tive massive galaxies in the (JFrei et aLlllQOq ) sample, while 
those of the Dwarf sample are of the order ~10^° M© (more 
akin to the Small Magellanic Cloud). We include one addi- 
tional dwarf (labelled DGILT) which was simulated using 
the same parent dark matter halo as one of the canonical 
dwarf simulations (DGl), but employed a lower star forma- 
tion threshold of 0.1 cm~^ (similar to the UW sample). 

We used SUNRISE |Jonsson 2006,), a Monte Carlo ra- 
diative transfer code, to generate a posteriori, mock images 
of the simulations, incorporating the effects of dust attenu- 
ation and spanning a range of viewing angles from edge- to 
face-on perspectives, as listed in Table [l] 



All galaxies are presented at inclinations of both 0° and 
90° , and in most cases also at a uniform distribution of inter- 
mediate inclinations. Within this distribution there are some 
correlations due to the intrinsic features of each galaxy that 
are maintained at all inclinations; e.g. in Figure 1111 DG2 
is clumpy at all inclinations. We view each galaxy from a 
range of inclinations, which allows us to compare them to 
observed galaxies, which are viewed at a range of inclina- 
tions. Furthermore, it increases the size of the sample given 
the limited number of galaxy simulations available at suf- 
ficient resolution. Subsequently we have multiple copies of 
the same physical system, thus the images are not truly in- 
dependent and care should be taken not to over-interpret 
any statistical conclusions. In Figures 1 to 4, SDSS r-band 
images of representative edge-on and face-on galaxies drawn 
from the observational (Figure 1) and simulated (Figures 
2-4) samples are shown. 



3 THE CAS STRUCTURAL PARAMETERS 

The three CAS p arame ters were recreated after 
IConsehce et al.l ll2000l) and IConseUcd (|2003l 'l. and cali- 
brated on the iFrei et al.l ([199y) sample of observed nearby 
galaxies in the 'r' and 'R' bands. Prior to the computation 
of each parameter, each galaxy was background subtracted 
by initially taking the median count value of a section at 
the edge of each galaxy's CCD frame, and subtracting this 
from each pixel within the array. The surface brightness 
profile was then inspected to ensure that the correct level 
of background was removed. 

Prior to computing the CAS indi ces, a characteri s tic ra - 
dius for each galaxy was determined. IConselice et al.l (|2000r ) 
recommends the use of the characteristic 77 radius which is 
1.5 rpctrosian- The Petrosian radius is defined as the radius 
r at which the ratio of the local surface brightness (in an 
annulus from 0.8r to 1.25r) to the mean surface brightness 
within r is 0.2. This is found by calculating the aforemen- 
tioned ratio for radii within an annulus 0.8 r to 1.25 r. 
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3.1 Concentration (C): 

The concentration parameter is a measure of how concen- 
trated (or difTuse) the central bulge component is with re- 
spect to the total flux of the galaxy. Concentration has 
a strong correlation with colour, and is related to other 
intrinsic features su ch as velocity dispersion and mass 
l|Graham et al.ll200ll ). The concentration of a galaxy is de- 
termined by locating the radii of the circular apertures that 
contain 20% (^0.2) and 80% (ro.s) of the total flux of the 
galaxy, and then taking their ratio: 



C = 51og 



ro.8 
ro.2 



(1) 



3.2 Asymmetry (A): 

The asymmetry parameter (A) is a measure of how sym- 
metric (or asymmetric) a galaxy is with respect to the total 
flux of the galaxy. Galactic asymmetries are predominantly 
caused by star formation, mainly in underlying structure 
such as spiral arms. Asymmetry is al so a strong indica tor of 
merger history, as demonstrated by IConselicd ([2003) . The 
asymmetry parameter initially involves locating the centre 
of a galaxy image, which is defined as the global minimum of 
asymmetry. The cent re was obtained using the centering al- 
gorithm specified bv lConselice et al.l ((20001), which involves 
an iterative trial-and-error process. An initial guess for the 
centre of the galaxy is made and the asymmetry for the 
specified pixel, and the surrounding eight pixels is deter- 
mined. The pixel corresponding to the lowest asymmetry is 
then redefined to be the centre. In order to locate the global 
asymmetry minimum, which corresponds to the centre, it is 
required that the process be repeated until the central pixel 
is also the pixel with the lowest value of asymmetry. 

The asymmetry itself is calculated by rotating the im- 
age 180° about the defined central point and subtracting 
the rotated image from the initial image, in order to de- 
termine the absolute value of the residuals. Following this, 
the residuals are divided by the quantity corresponding to 
the summation of all the pixels in the galaxy. Subsequently, 
a similar method is then applied to a small section, taken 
from a corner of the image which contains background only. 
This quantity is then multiplied by a scale factor so that it 
is comparable to the size of the galaxy image. This quantity 
is divided by the summation of the total flux of the original 
galaxy image. This is done in order to determine the asym- 
metry of the background so that it may be discounted from 
the galactic asymmetry. The calculation for the asymmetry 
is as follows: 



A = 



E ^ - -fiso T,B-Bi 



EI 



EI 



(2) 



where / is the image flux, hso is the rotated image flux, B 
is the background, and Biso is the rotated background. 



3.3 dumpiness (S): 

The measure of high-spatial frequency dumpiness (S) is di- 
rectly related to a galaxy's star formation and high fre- 
quency spatial power. The dumpiness is deflned as the 
patchiness of the high fre quency light distribution within 
a galaxy (|Conselicel 120031 ). It is computed by subtracting 



a degraded-resolution version of the galaxy image from the 
original galaxy image. As with the asymmetry parameter 
the method is repeated on a section of sky background in 
order to isolate the fraction of the signal that is due to the 
galaxy. Again, a corner of the image containing background 
only is selected, replicated, and smoothed, so that an iden- 
tical procedure may be carried out on both the background 
and the galaxy image. 

As the central bulge component is a prominent feature 
in the simulated galaxies, it contributes to a large portion 
of the high frequency structure and hence dumpiness of 
the simulated galaxies. As it is already a known issue that 
simulated galaxies tend to have large central bulge compo- 
nents, with respect to their observed counterparts, it was 
decided that the standard central region specifled by Con- 
selice (2003) would be extended to include the total central 
bulge component of each galaxy. Furthermore, the inferred 
total flux of the galaxies would also exclude the central re- 
gion, in order to avoid any bias due to the varying sizes of 
the bulge components of the galaxies. The size of the cen- 
tral bulge component was determined from the light proflle 
of the galaxy by considering the radius at which the proflle 
changed from ln(/) — r^'*, to exponential decay. The im- 
age was then examined to ensure that the correct area had 
been removed, prior to the computation of the dumpiness 
parameter. 

The calculation for the dumpiness is as follows: 



5 = 10 X 



EI -Is EB-Bs 



EI 



EI 



(3) 



where / is the image flux. Is is the smoothed image flux, B is 
the background, and B3 is the smoothed background. Fur- 
ther details on the derivation of the dumpiness parameter 
are provided in Appendix A. 



4 COMPARISON OF OBSERVATIONS WITH 
SIMULATIONS 

The CAS parameters form a well-established system that 
has been used on occasion to analyse galaxy simulations, 
although not i n any systematic manner. For example, 
(|Conselicell2006l ) applied th e system to dark matter merger- 
remnant simulations, and (JLotz et al.ll2008h extended this 
to the analysis of equal-mass gas-rich disc merger remnants. 
Our unique contribution to this area is, for the first time, an 
application to a suite of simulations realised within a fully 
cosmological framework. Such a comparison should enable 
underlying physical differences to be revealed and allow for 
scaling relations to be tested using simulations in a new and 
unique manner. The main purpose of this paper is to con- 
sider the morphological details of the galaxies produced in 
simulations, and uncover new information regarding the dif- 
ferences between observed and simulated galaxies. To avoid 
simply re-identifying known issues and biases, we have mod- 
ified slightly the underlying CAS parameters, as explained 



4.1 Concentration (C): 

The overproduction of the centrally-concentrated bulge 
component of most simulated disc galaxies (e.g., Stinson 
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Table 1. CAS parameters for simulated galaxies, including values 
obtained at different inclinations. The names of the galaxies were 
defined by the teams that simulated them (see text for details); 
the specific names are arbitrary for the purposes of this study. 



UW Sample 


Galaxy Inclination (°) 


C(r) 


A(r) 


S{r) 


Gall 





3.49 


0.351 


0.485 




90 


3.26 


0.312 


0.680 


h277 





4.96 


0.215 


0.126 




90 


3.82 


0.161 


0.508 


MW1024 





3.97 


0.180 


0.181 




90 


3.65 


0.180 


0.357 


Dwarf Sample 


Galaxy Inclination (°) 


C(r) 


A(r) 


S{r) 


DGl 





3.38 


0.385 


0.378 




45 


3.38 


0.232 


0.453 




90 


3.47 


0.369 


0.667 


DGILT 





2.80 


0.267 


0.259 




30 


3.01 


0.281 


0.253 




45 


2.87 


0.275 


0.303 




60 


3.18 


0.260 


0.380 




90 


4.06 


0.231 


0.761 


DG2 





3.66 


0.574 


0.217 




30 


3.17 


0.570 


0.193 




45 


3.19 


0.568 


0.228 




60 


3.18 


0.552 


0.197 




90 


3.17 


0.376 


0.699 


DG3 





3.28 


0.315 


0.139 




30 


3.31 


0.317 


0.126 




45 


3.29 


0.314 


0.134 




60 


3.22 


0.311 


0.177 




90 


3.48 


0.200 


0.664 


MUGS Sample 


Galaxy Inclination (°) 


C(r) 


A(r) 


S(r) 


gl536 





4.51 


0.158 


0.078 




45 


4.46 


0.174 


0.084 




90 


4.11 


0.123 


0.201 


g21647 





3.49 


0.162 


0.104 




45 


3.31 


0.150 


0.061 




90 


3.49 


0.173 


0.314 


g22437 





3.71 


0.240 


0.258 




30 


3.49 


0.193 


0.186 




45 


3.34 


0.209 


0.159 




60 


3.31 


0.222 


0.146 




90 


3.45 


0.209 


0.192 


g24334 





3.57 


0.235 


0.362 




30 


3.44 


0.235 


0.424 




45 


3.44 


0.230 


0.406 




60 


3.36 


0.215 


0.393 




90 


3.54 


0.246 


0.457 


g25271 





3.63 


0.106 


0.007 




45 


3.63 


0.086 


0.007 




90 


4.06 


0.148 


0.100 


g3021 





3.71 


0.206 


0.277 




30 


3.76 


0.216 


0.324 




45 


3.84 


0.235 


0.393 




60 


3.91 


0.234 


0.407 




90 


4.03 


0.209 


0.419 


g5664 





3.16 


0.146 


0.359 




45 


2.93 


0.170 


0.330 




90 


3.39 


0.142 


0.726 


g7124 





4.51 


0.077 


0.000 




/li; 


A nn 


n noQ 


n noQ 



et al. 2010, and references therein; c.f., Governato et al. 
2010), is a long-standing problem associated with cosmo- 
logical simulations of galaxy formation. By definition, the 
concentration parameter is a proxy for the dominance of the 
central region, simulation or not, although consideration of 
the surface brightness profiles suggested that the basic def- 
inition not be altered here. Instead, we concentrate here on 
the relationship between the size of the central component 
and disk of each galaxy, and compare this with observations. 
Another discrepancy between observed and simulated 
galaxies is their noise content as noise does not overlay the 
simulated galaxies. Unlike the other CAS parameters, the 
measure of concentration does not include a background sub- 
traction stage. This is due to the negligible effects of noise 
on the measure of concentration. For this reason, no adjust- 
ments were made to the concentration parameter regarding 
background noise. 

4.2 Asymmetry (A): 

The measure of asymmetry considers the residuals produced 
when a galaxy image is subtracted from the same image 
that has been rotated by 180°. The central bulge compo- 
nent of a galaxy is generally symmetric, and for this reason 
no adaptation was required to account for the excessive cen- 
tral bulge component of the simulated galaxies. However, 
it was discovered that the effects of noise were critical, as 
the background component was comparable in magnitude 
to the component due to the galaxy itself. It was imperative 
to ensure that this aspect of the algorithm was correct, as 
simulated galaxies contain essentially zero background noise 
compared to their observed counterparts, the latter of which 
contain varying amounts of noise depending upon exposure 
time, sky background, and distance. For this reason, any dis- 
crepancies concerning the computation of the background 
would lead to an inaccurate comparison. 

In order to test this hypothesis, an artificial array of 
noise was created, and the galaxies were analysed both be- 
fore and after the application of this noise array. We repli- 
cate the signal-to-noise values of the Frei et al. (1996) sam- 
ple by first pixelating the simulated images to the same 
resolution, then scaling the flux to be comparable to that 
in the images of the real galaxies, and finally adding a 
normally-distributed noise array of sufficient amplitude to 
reproduce the signal-to-nois e of the observed galaxies in Frei 
at al. (1996). Furthermore, IConselice et al.l (|2000h demon- 
strated that accurate asymmetry calculations can only be 
performed, without corrections, on galaxies with signal-to- 
noise greater than 100. For this reason, each galaxy was 
tested individually to ensure that the application of these 
characteristic noise values led to signal-to-noise ratios in ex- 
cess of 100. 

A comparison of the asymmetry results derived from 
the simulated galaxies, both before and after the applica- 
tion of noise, can be seen in Figure [5] It was found that the 
inferred asymmetry values differed only marginally after the 
application of noise (in large part because of the stringent 
signal-to-noise limit imposed), although the impact varied 
for each galaxy and appeared to be more substantial for 
more asymmetric galaxies. Whilst there is some variation 
visible in Figure [5] due to the application of noise, it is not 
significant. However, to ensure that a fair comparison was 
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Figure 5. Relationship between the asymmetry of simulated 
galaxies both before and after the application of noise. The solid 
line represents the least squares regression line for which the equa- 
tion is j; = 1.15x — 0.002, and the dashed line corresponds to the 
one-to-one relation. Here, the points are coloured to represent 
the different galaxy samples: MUGS = purple, UW = green and 
Dwarf = red. The choice of symbol corresponds to galaxy inclina- 
tions of to 45 degrees ('face-on': diamonds) and 45 to 90 degrees 
('edge-on': asterisks). 



made between observed and simulated galaxies, it was de- 
cided that noise would be applied to the simulated galaxies 
prior to the computation of the asymmetry parameter. 



4.3 dumpiness (S): 

When initially computed, the central bulge component dom- 
inated the dumpiness calculation for the simulated galax- 
ies, suggesting that the simulations were excessively clumpy 
with respect to their ob served counterp arts. Prior to cal- 
culating the dumpiness, ()Conselicell2003l ) forces the central 
region, which corresponds to 0.05 times the Petrosian radius, 
to zero, to ensure that any unresolved region of the central 
bulge component is excluded. For the purpose of this study 
this region has been extended to encompass the total cen- 
tral bulge component of each galaxy, thus allowing for the 
comparison of observed and simulated galaxies outside the 
bulge region, which we know is overproduced in these sim- 
ulations. By making this alteration, the dumpiness of the 
disk region only is calculated, excluding any high frequency 
structure associated with the central bulge component. 

As the area of the central bulge differs for each galaxy, in 
order to avoid any bias due to bulge size, it was then decided 
that the central bulge component would also be excluded 
from the summation of the total flux which is present in 
the denominator of the equation for the computation of the 
dumpiness (recall, i)3.3p . 

The application of noise had a negligible effect upon the 
values generated by the dumpiness algorithm. This suggests 
that the subtraction algorithm for the dumpiness parameter 
is more than adequate for minimising the impact of noise on 
its calculation. For this reason it was deemed unnecessary 
to apply noise to the simulated galaxies in the computation 
of the dumpiness. 



5 DISCUSSION 

For each simulation, surface brightness profiles, for a range 
of inclinations, were generated. The simulations were 'pix- 
elated' to a scale corresponding to 100 pc/pixel, which, as 
shown by Consclicc ct al. |200Ci), appears adequate for CAS 
analysis. The iFrei et al.l (|l996r ) sample of observed, nearby 
galaxies, used in this study, all possess resolutions compa- 
rable to (or better) than the same 100 pc/pixel, and hence 
are again ideal for use with the CAS system. As all the 
galaxies used in this study are above the threshold value of 
~1 kpc/pixel, wh ere the accuracy of th e CAS parameters 
begins to decline (|Conselice et al.ll2000f ), a comparison be- 
tween simulated and observed galaxies can be made with a 
high deg ree of confidence. 

The iFrei et al.l l|l996l) sample of galaxies contains only 
nearby high surface brightness galaxies. It is not meant to 
be a complete and/or unbiased compilation, as it mostly 
excludes irregular and low surface brightness dwarf galax- 
ies, which contribute significantly to the population of lo- 
cal galaxies. Similarly, the sample of simulated galaxies is 
incomplete in the sense that it does not contain any irreg- 
ular or elliptical galaxies, resulting in a restricted range of 
colours; this impacts particularly on the re d end of the distri- 
bution, where we lack counterparts to the lFrei et al.l (|l996r ) 
sample. This does place constraints on the ability to com- 
pare observed and simulated galaxies, and is considered dur- 
ing this analysis by incorporating Hubble type in many of 
the figures. 



5.1 Concentration 

In Figure |6] both observed and simulated galaxies demon- 
strate a strong correlation between colour and concentra- 
tion. This trend naturally arises for observed galaxies where 
early-type galaxies are redder in colour and also tend to be 
more centrally-concentrated, whereas later-types are bluer 
due to ongoing star formation and tend to be less centrally- 
concentrated. These values are well quantified by Hubble 
type and hence, for the observed galaxies, it can be seen 
that the different symbols that represent the various Hubble 
types are contained within well-defined regions. The visual 
appearance of a trend among the simulated galaxies comes 
more from the systematic offset between the Dwarf sample 
in red and the MUGS and UW samples in purple and green. 
While there is a trend within the MUGS sample of galaxies, 
it is not nearly as visually striking. 

It is evident from Figure]^ that the simulated galaxies 
tend to be bluer than the observed galaxies, both in general 
and with respect to their Hubble type. The main cause of 
this blueward offset is the presence of excessive late-time 
star formation in the simulated discs. Furthermore, there 
are a lack of elliptical galaxies in the simulated sample, with 
the earliest galaxies being lenticular. However, in a relative 
sense, there is a definite relationship between colour and 
Hubble type, as expected. 

The simulated galaxies in Figure |S] show a similar range 
of concentrations to those observed for real galaxies, and also 
appear to be divided into groups by Hubble type. Further- 
more, the simulated galaxies demonstrate a positive linear 
trend that, although offset due to colour differences, is com- 
parable to the trend demonstrated by the observed galaxies. 



Advanced Morphological Galaxy Classification 7 



n E - 


SO 


X So 


- Sb 


O Sc 


- Sd 



0.2 1 

2.5 



'* -M 

""^ t 



}>\, 



a 



a 



n ■» 






3.5 4.0 

Concentration 



1 D 
S < 



0.2 1 

0.0 



*; □ 



>n 



W. 



0>* : 



^■ 






E - 


S'",' D 


So - 


- Sb X 


Sc - 


- Sd O 



Figure 6. Relationship between concentration and colour (B-V) 
for both observed and simulated galaxies. The points are coloured 
to represent the different galaxy samples: observed = blue (and 
slightly larger in size relative to the other symbols), MUGS = 
purple, UW = green, and Dwarf = red. The symbols correspond 
to the classical Hubble types, as noted in the inset. 



Concentration is also related to surface brightness; it is ex- 
pected that the low concentration values generated by the 
majority of the simulated dwarf galaxies is primarily due to 
their low surface brightness. Furthermore, the dwarf galax- 
ies in this sample are all younger disc-like galaxies which are 
observed to be generally more diffuse. Other factors which 
influence the concentration parameter include the mass and 
velocity dispersion, both of which are low for the dwarf 
galaxy sample. 

The MUGS sample of simulated galaxies in Figure[6]ap- 
pears to extend to higher values of concentration than both 
the UW sample of Milky Way-like galaxies and the Dwarf 
sample. This is most likely due to the presence of earlier-type 
discs that have greater masses within the unbiased MUGS 
sample (c.f., with the UW sample, which were chosen a pri- 
ori to more closely resemble later- type discs). Moreover, the 
MUGS sample of galaxies have a higher density threshold 
for star formation (1 cm"''), which could lead to preferential 
star formation in the denser (central) regions of the galaxy 
(Agertz, Teysier & Moore 2010), although this behaviour is 
not seen in dwarf galaxy simulations (Governato et al. 2010). 
The MUGS sample also has a greater average circular ve- 
locity (150 - 250 km/s), compared to the other simulated 
samples of galaxies, whi ch is a parameter k nown to corre- 
late with concentration l|Graham et al.ll200ll ). 



5.2 Asymmetry 

In Figure [71 the observed galaxy sample demonstrates a 
strong anti-correlation between asymmetry and colour. This 
correlation stems from the bluer, younger galaxies gener- 
ally containing spiral arms with dust lanes and star form- 
ing regions, and the older, elliptical galaxies generally being 
devoid of structure. The simulated galaxies follow an ap- 
proximately identical trend to that of the observed galaxies, 
however, due to the bluer nature of the selected simulations, 
they form a blu e extension of the observed galaxy trend. 
IConselicd ((2003) showed that there was a distinct region in 



Figure 7. Relationship between asymmetry and colour (B-V) for 
both the observed and simulated galaxies. Symbols and colours 
are the same as in Figure |6] 
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Figure 8. Relationship between the asymmetry (A) of simulated 
galaxies, without the application of noise, and colour (B-V). Sym- 
bols and colours as in Figure |5] 



the bottom left corner of the colour-asymmetry plot where 
no galaxies were bluer and more symmetric. This region is 
well-defined by the simulations which demonstrate a very 
strong correlation between colour and asymmetry. 

The earlier-types within the MUGS sample are redder 
in colour with respect to the other simulations, and tend 
to be more symmetric. Again there is a strong distinction 
between the various Hubble types, which is similar to the 
correlation demonstrated by the observed galaxy sample. 
However, due to the lack of redder simulated galaxies with 
low, approximately zero, asymmetry values, this trend does 
not extend to the upper left corner of Figure [T] where the 
redder, symmetric, galaxies reside in nature. Again, this can 
be traced to the ongoing late-time star formation in the 
earlier-type simulations which is not seen in nature. In Fig- 
ure [7] some of the dwarf simulated galaxies appear to be ex- 
cessively asymmetric relative to the other simulations. The 
large variation of apparent flux levels within these systems 
is most likely the cause of this increased in asymmetry. 

IConselice et al.l (|200G| ) showed that it is almost always 
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Figure 9. Relationship between dumpiness (S) and colour (B- 
V) for both the observed and simulated galaxies. Symbols and 
colours as in Figure |6] 



the MUGS sample that possess approximately zero dumpi- 
ness, consistent with their observed counterparts. 

In Figure]^ it can be seen that there is a group of galax- 
ies (mainly from the Dwarf sample, with the exception of one 
MUGS and one UW galax;y), that are excessively clumpy. 
All the simulated galaxies that are excessively clumpy also 
have inclination angles of 90° (see Table [l]). Thus a study 
into the effects of dust attenuation on the dumpiness pa- 
rameter was carried out. It was found that the dumpiness 
of the simulated galaxies was not affected by the contri- 
bution of the dust from SUNRISE (|jonssonll200^ . It was 
therefore concluded that the Dwarf sample, and Gall from 
the UW sample, are simply intrinsically clumpy. In the case 
of the Dwarf sample, this is consistent with an analysis of 
the power spectru m of the cold gas ass ociated with their 
interstellar media (IPilkington et al.l 120111 ). This may be at- 
tributed to the relatively low surface brightnesses associated 
with the dwarfs, as this causes the high frequency residuals 
to dominate in the calculation of the dumpiness. 



the inclined galaxies that contribute to the scatter in the 
colour-asymmetry diagram and that the less inclined galax- 
ies are systematically higher. In Figure [S] the edge-on and 
face-on galaxies have different symbols, and from this it can 
be concluded that while the inclined, observed galaxies do 
tend to be somewhat more significant outliers, in detail, the 
scatter of the simulated galaxies is not a function of inclina- 
tion angle. Also, the MUGS and UW galaxies that are less 
inclined do appear to be systematically higher in the plot, 
as found for the observed galaxies. 

Noise was applied to the simulated galaxies to ensure 
that a fair comparison was made when considering the asym- 
metry. The addition of a noise array to the simulated galaxy 
images had the benefit of allowing for a more accurate com- 
parison with the observed galaxies. In Figure [8] the range 
of asymmetry values, without the addition of noise, can be 
seen. The absence of noise only appears to affect the out- 
come of the asymmetry parameter at greater asymmetry 
values by restricting the range of asymmetries. 



5.3 dumpiness 

In Figure|9l it can be seen that there is an inverse correlation 
between colour and dumpiness for both the observed galaxy 
sample and the MUGS simulated galaxy sample, although 
the anti-correlation of the MUGS galaxy sample is much 
weaker (and shifted in colour). The UW sample shows a 
positive correlation (though over a limited colour range) and 
the Dwarf sample shows no correlation. As the earlier-type 
disc galaxies are generally devoid of star forming regions 
and internal structure, and also tend to be red in colour, 
it follows that they would generally be located at the top 
left of Figure |9l where the smoother, redder galaxies reside. 
Conversely, later-type galaxies tend to be bluer and contain 
high frequency residuals. For this reason, the anti-correlation 
is the expected correlation for Figure [9l 

The simulated and observed galaxies demonstrate a 
strong correlation between dumpiness and Hubble type. Due 
to the lack of redder galaxies, as with asymmetry, the sim- 
ulated galaxies do not reach the top left corner of Figure (9] 
However, unlike asymmetry, there are several galaxies from 



5.4 CAS comparison 

The CAS parameters are based solely on the morphology of 
galaxies, but are directly related to many intrinsic qualities. 
By comparing CAS parameters amongst themselves, as op- 
posed to CAS versus colour, sampling biases, such as the 
colours of stellar populations, may be avoided. This allows 
for the structures of the galeixies to be compared directly, 
which may further emphasise the differences between ob- 
served and simulated galaxies. 

In Figure [TD] the observed galaxies from the lFrei et al.l 

(|l996r ) sample demonstrate a general negative correlation 
between concentration and dumpiness. Galaxies that are 
smooth tend to be concentrated and those that are clumpy 
tend to be more diffuse. This correlation is generated be- 
cause younger, disc galaxies tend to be both clumpy due 
to active star formation and diffuse, while older elliptical 
galaxies tend to be smoother and more dense. 

The majority of simulated galaxies lie within the nor- 
mal range of observed galaxies, in Figure ITOl It is apparent 
though that some of the simulations are excessively clumpy 
in comparison to their observed counterparts. As can be gar- 
nered from Table 1, all these simulated galaxies which are 
clumpy outliers only occur when viewed at an inclination 
of 90°. This group of galaxies display approximately aver- 
age concentration values, whilst being excessively clumpy, 
in comparison to both their observed and simulated galaxy 
counterparts. In regards to the excessively clumpy dwarfs, 
we note that the high frequency structure present in these 
simulated galaxies is consistent with the power spectrum 
analysis of the associated ISM by Pilkington et al. (2011). 

In Figure 1101 the MUGS sample demonstrates an ap- 
proximately equivalent trend to the observed galaxy sample, 
despite their excessive bulge sizes. This is unexpected, as the 
central bulge component is disregarded in the computation 
of the dumpiness, but is included for the concentration pa- 
rameter. However, with respect to the Dwarf and UW sam- 
ples, the MUGS sample contains galaxies that are generally 
more concentrated. It is likely that the greater mass, veloc- 
ity dispersion, and higher star formation threshold of the 
MUGS sample all contribute to both the larger bulge sizes 
and concentration of these galaxies, with respect to the other 
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Figure 10. Relationship between tlie concentration (C) and 
dumpiness (S) parameters for both the observed and simulated 
galaxies. The points are coloure d to represent the different galaxy 
samples. The iFrei et al.l l|l996h sample of observed galaxies are 
coloured blue with larger symbols relative to the simulated sam- 
ples. The MUGS sample are shown in purple, the UW sample in 
green, and the Dwarfs represented with red symbols. The symbols 
denote the inclination angle, as in Figure [5] 
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Figure 11. Relationship between dumpiness (S) and asymmetry 
(A) parameters for both observed and simulated galaxies. Sym- 
bols and colours as in Figure [TO] 



simulated galaxy samples. Furthermore, it is likely that these 
characteristics produced the earlier-type discs in the MUGS 
sample, which are expected to be less clumpy in nature. 

In Figure llll it can be seen that there is a strong pos- 
itive correlation between the dumpiness and asymmetry of 
observed galaxies. This correlation arises from the later-type 
galaxies tending to be clumpy and asymmetric, due to their 
star forming regions and intrinsic structure, and earlier-type 
galaxies tending to be devoid of structure and star form- 
ing regions, which pertains to them being more smooth and 
symmetric. The simulated galaxies appear to follow the same 
correlation. However, the sample of highly symmetric ob- 
served galaxies are not reproduced within the simulations. 

In Figure llll the Dwarf sample of galaxies form three 
groups; those which lie within the normal range, those which 
are excessively clumpy, and those that are excessively asym- 



Figure 12. Relationship between concentration (C) and asym- 
metry (A) parameters for both observed and simulated galaxies. 
The points are coloured to represent the different galaxy samples. 
Symbols and colours as in Figure [Tol 



metric. As previously discussed, those galaxies that are ex- 
cessively clumpy all have 90° inclination angles. Figure [TT] 
shows that the high frequency spatial dumpiness of these 
galaxies is not directly related to asymmetry, as those galax- 
ies that demonstrate excessive dumpiness possess a wide 
range of asymmetries with respect to other galaxies from 
the same sample. 

The group of simulations which are seemingly highly 
asymmetric outliers at the top of Figure [TT] is in fact just 
one simulated galaxy (DG2) viewed at various inclinations. 
This particular galaxy has a large star forming complex at 
z = which is offset from the centre of the galaxy (see 
Governato et al. 2010), explaining both the dumpiness and 
asymmetry of this galaxy. 

In Figure [T2l there is an apparent negative correlation 
between asymmetry and concentration for both observed 
and simulated galaxies. This can be attributed to the fact 
that later-type galaxies are both diffuse and asymmetric, 
and earlier-type galaxies are generally dense with less inter- 
nal structure. Within the Dwarf galaxy sample there is again 
a single galaxy viewed at various inclinations that demon- 
strates low concentration values coupled with extremely high 
asymmetry values. The large star forming complex in this 
simulation has not resulted in a concentration being outside 
the region defined by observed galaxies. 



6 CONCLUSION 

We have applied measures of concentration (C), asymmetry 
(A) and dumpiness (S) to a sample of simulated galaxies 
which have a range of masses and morphologies, as well as 
to a sample of observed galaxies. We have explored the cor- 
relations between these parameters as well as between each 
of them with B-V colour. In general, reasonable agreement 
was found between the simulated and observed populations 
in these relations, although some differences become appar- 
ent, as summarised below. 

Overall, the trend generated by the observed galaxy 
sample, with respect to colour (B-V), was approximately 
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replicated by the simulated galaxies, although there is a 
distinct lack of redder simulated galaxies, which is re- 
lated to a well-established inability of these simulations 
to halt late time star formation in early type galaxies 
l|Kawata fc GibsonllJOOsh . 

It was established that Hubble type is strongly linked 
with the three CAS parameters for the both the observed 
and simulated galaxy samples. This is expected as later- 
type galaxies, such as Hubble type Sc-Sd, have star form- 
ing regions, resulting in significant internal spatial struc- 
ture, which manifests itself in significant asymmetries and 
dumpiness, and possess low concentration values. Con- 
versely, earlier- type galaxies are smoother, more symmet- 
ric, and more concentrated, due to their relative lack of star 
forming regions, structure, and their evolutionary histories. 

The concentration parameter was found to be robust 
and easily replicated. The values obtained for the simu- 
lated galaxies demonstrated an approximately equal range 
of values compared with their observed counterparts. The 
fewer number of simulated galaxies with high concentration 
merely refiected the relative paucity of early-type galaxies 
compared with the observed sample. 

It was found that the range of asymmetries of the sim- 
ulated galaxies did not extend to as low values as those of 
the most symmetric observed galaxies. This may be related 
to the ongoing star formation in the simulated galaxies, the 
same process which results in a dearth of red galaxies in the 
simulations. It is generally the early-type observed galax- 
ies which are highly symmetric. Nevertheless, it remains a 
challenge for simulators to create highly-symmetric galaxies 
within a hierarchical structure formation paradigm. 

A single dwarf galaxy simulation also shows excessive 
asymmetry compared to the observed sample. This is due 
to the existence of a large star forming complex offset from 
the centre of that simulation. Detailed exploration of the 
resolution of star forming complexes and star clusters within 
simulations is beyond the scope of this study, but it appears 
that the physical size of the star forming complex in at least 
one of the simulations is larger than observed, at least within 
the Frei a et al. (1996) sample. 

As it is a well known issue that simulated galaxies have 
dispr oportionately large bulge components IjStinson et al.l 
[201y), it was decided that a central region, pertaining to 
the size of the bulge component, would be removed for all 
galaxies in the calculation of the dumpiness. Following this, 
we found that the simulated and observed galaxies showed 
similar spreads in the values of dumpiness. The early-type 
disc galaxies from the MUGS sample, which have lower 
rates star formation occurring in the disk at 2 = 0, have 
low S, while the later- type simulated galaxies continue to 
have high levels of star formation, and high values of S, 
consistent with their observed counterparts. Analysis of the 
dumpiness values generated by the simulated galaxies also 
demonstrated that several simulated galaxies were exces- 
sively clumpy when viewed at an inclination angle of 90°. 
At least in relation to the dwarfs, this is likely related to 
high frequency power of these galaxies, when analysed at 
an edge-on inclination, a characteristic that is not consis- 
tent with observations (cf., Pilkington et al. 2011, for the 
high-frequency analysis of the ISM of these dwarfs). 
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APPENDIX A: COMPUTATION OF THE HIGH 
SPATIAL FREQUENCY CLUMPINESS 
PARAMETER (S) 

For the purpose of comparing simulated galaxies with those 
observed, the central bulge component of each galaxy was 
removed in the computation of the high frequen cy spatial 
dump iness. However, the method proposed by IConselicd 
l|2003l ) favoured the removal of a central region, amounting 
to 0.05 X the Petrosian radius (rpetrosian), for the purpose of 
excluding any unresolved regions of the image only. As the 
measure of high spatial frequency dumpiness specified by 
IConselicd (12003) forms an integral part of the CAS system, 
the method of computation is presented here. 

Initially the image is sky-subtracted, flat-fielded, and 
the characteristic radius determined. The radius is deflned 
as 1.5 X rpctrosian (further details can be found in ?). If the 
data is undersampled, the original image undergoes a box- 
car smooth of 5x5 pixels prior to any further computation; 
this image is then re-deflned as the original image. The next 
step involves smoothing the image by a factor r-pctrosian/6, 
again using a boxcar smooth. Following this, a map of the 
residuals is created by subtracting the smoothed image from 
the original image. Within this residual image, all the nega- 
tive pixels, and the pixels within the central region of 0.05 x 
T'potrosian, are then forced to zero, leaving behind the high fre- 
quency residuals only. The total flux from the residual map, 
within the specifled radius, is then summed and divided by 
the total flux of the initial galaxy, within the specified ra- 
dius. The value generated equates to the dumpiness of the 
galaxy including the dumpiness in the background. 

To ensure that the results of the dumpiness parameter 
determination are unaffected by noise, further background 
correction is then required. This involves selecting an area 
from the image that contains background only and repeat- 
ing the aforementioned procedure. Note that the background 
should also be smoothed by a boxcar smooth of 5x5, if 
the data is undersampled. Once selected, the area should 
be smoothed by a boxcar smooth of rpetrosian/6, using the 
value of the Petrosian radius obtained from the galaxy im- 
age. A background residual map is then created by sub- 
tracting the smoothed background image from the original 
background image and forcing the negative pixels to zero. 
The background flux is then be summed and normalised to 
the the number of pixels that constitute the total galaxy, 
within the specified radius. This value is then divided by 
the summation of the total flux of the galaxy thus produc- 
ing the dumpiness in the background. The dumpiness in 
the background is then subtracted from the aforementioned 
dumpiness in the galaxy, which is computationally defined 



